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Dissociation of O2 at Al(lll): The Role of Spin Selection Rules 
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A most basic and puzzling enigma in surface science is the description of the dissociative adsorption 
of O2 at the (111) surface of Al. Already for the sticking curve alone, the disagreement between 
experiment and results of state-of-the-art first-principles calculations can hardly be more dramatic. 
In this paper we show that this is caused by hitherto unaccounted spin selection rules, which give 
rise to a highly non-adiabatic behavior in the 02/Al(lll) interaction. We also discuss problems 
caused by the insufficient accuracy of present-day exchange-correlation functionals. 

PACS numbers: 82.20.Kh, 82.20.Gk, 68.35.Ja 
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Oxygen-metal interactions are responsible for every- 
day phenomena like corrosion, and form the atomic-scale 
basis behind numerous technological applications like ox- 
idation catalysis. It is therefore most discomforting that 
despite several decades of research in surface science, the 
initial step in the oxygen-metal interaction, namely the 
dissociation process of O2 molecules over metal surfaces, 
is not yet understood. This is in particular so for what 
is often called the most simple metal surface, namely 
Al(lll): a close-packed surface of a nearly-free electron 
metal. For the initial interaction of Oawith Al(lll) ex- 
periments have consistently shown P, Q that the initial 
dissociative sticking probability for thermal O2 is very 
low (about 2%). Density- functional theory (DFT) cal- 
culations, on the other hand, found that dissociation is 
not hindered by energy barriers which implies that 
the initial sticking coefficient should be very high (about 
100%). Another intriguing aspect of the 02/Al(lll) sys- 
tem is that at very low coverages the distribution of ad- 
sorbed oxygen atoms is random, even when adsorption 
is performed at temperatures at which thermal diffusion 
can not play a significant role Thus, it is impossible 
to trace back which two adatoms stem from the same 
molecule. Initially this led to the suggestion that the ad- 
sorption energy is used to trigger the diffusion of "hot 
adatoms" 01 ■ More recently, a different explanation has 
been suggested ("abstraction"), where only one 0-atom 
is adsorbed and the other one is repelled back into the 
vacuum Q. Again, theoretical work, so far, does not give 
a clue why this may be so. Thus, one may ask, what we 
can trust in surface science when understanding of such 
a most basic and simple system for molecule-surface in- 
teractions is so clearly lacking. 

Figure n summarizes the experimental data for the ini- 
tial sticking coefficient as function of the kinetic energy 
of incoming O2 molecules for a molecular beam at nor- 
mal incidence (full diamonds) I'l, as well as the result, 
of what has hitherto been the standard theoretical treat- 
ment (labeled as "theory-adiabatic" ) . Also shown is the 
result of the approach taken in the present paper (labeled 
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FIG. 1: Initial sticking curve of O2 at Al(lll), based on 
the adiabatic (empty triangles) and the spin-triplet (empty 
circles) potential-energy surfaces using the RPBE functional. 
The experimental data (solid diamonds) are from ref. Q. 



as "theory-triplet"), which will be detailed below. Obvi- 
ously, there is hardly any similarity between the "theory- 
adiabatic" curve and the experimental result. Though 
we called this the "standard theoretical treatment", we 
note that already the calculations behind the "theory- 
adiabatic" curve (and also behind the "theory-triplet" 
curve) are much more elaborate and advanced than typ- 
ical approaches to obtain the initial sticking coefficient: 
All theoretical results presented in this paper were ob- 
tained from extensive all-electron DFT calculations using 
the DMoP code 0. This provided the six-dimensional 
potential-energy surface (PES) for the 02/Al(lll) sys- 
tem at more than 1500 geometries of the two oxygen 
atoms, keeping the substrate frozen. These PES data 
points were then interpolated by a neural- network Q , 
enabling us to perform molecular dynamics (MD) calcu- 
lations for about 100,000 trajectories, including all possi- 
ble initial molecular orientations. Thus, this approach Q 
grants a controlled and good statistics, in contrast to "on- 
thc-fly ah initio MD" , which gives (for a frozen substrate) 
the same trajectories, but where due to the high CPU 
cost at best only ~ 50 trajectories could be performed 
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even on todays biggest computers. 

Still, "on-the-fly ab initio MD" has the advantage that 
it can also be used beyond the frozen substrate approxi- 
mation. To check on the validity of our treatment, we 
therefore performed 24 ab initio MD runs, where the 
full dynamics of the Al surface atoms was taken into ac- 
count. These studies show that the adsorption energy 
is efficiently transferred to strong surface vibrations, and 
that the oxygen adatoms do not move far. Thus, the 
"hot adatom" concept is not supported. In all studied 
trajectories the Al(lll) surface got only affected, when 
the O2 was quite close to the surface, i.e. when 0-Al 
bonds were being formed and the 0-0 bond notably 
weakened (at molecule-surface distances below « 2.5 A). 
Before this point, the O2 trajectories were not changed 
by the substrate vibrations, and in particular all incom- 
ing O2 molecules are found to dissociate, fully confirm- 
ing the adiabatic result shown in Fig. ^ We also per- 
formed a systematic comparison using different exchange- 
correlation (xc) functionals, including the PBE Jlj and 
RPBE [13. The resulting PESs look different iiTsome 
details, however, the resulting sticking curve is always 
essentially the same as the "theory-adiabatic" curve in 
Fig. n Hence, neither the approximate xc treatment, 
nor the frozen substrate approximation can account for 
the dramatic disagreement between the theoretical and 
experimental results. We therefore conclude that the ori- 
gin must be more fundamental, namely in the assumed 
adiabatic description, restricting the impinging molecule 
to the electronic ground state of the combined O2 / Al sys- 
tem at each point of the O2 trajectory. Based on less rig- 
orous studies, this had been suggested previously 

Inspecting the six-dimensional adiabatic PES reveals 
immediately an obvious flaw of the adiabatic descrip- 
tion, independent of the employed xc functional: Even 
at largest distances the electron chemical potentials of 
the O2 molecule and the Al(lll) surface align, which 
is achieved by some electron transfer towards the O2 
molecule. Obviously, in reality charge transfer will occur 
only when the two systems are getting close for a suffi- 
ciently long period of time. Recently, Hellman et al. 
considered the influence of charge transfer by employ- 
ing an approach, where they replaced the Al(lll) sur- 
face by jellium and treated the kinetic-energy operator 
in the Thomas-Fermi- Weizsacker approximation. Then, 
two one-dimensional diabatic PESs were constructed, one 
where the O2 molecule was kept neutral and one where a 
full electron was transferred This description could 
indeed account for the qualitative shape of the experi- 
mental sticking curve, as could Binetti et al. (12], follow- 
ing a comparable approach, but considering four different 
diabatic model PESs. Both treatments point therefore at 
the possible importance of non-adiabatic effects, but due 
to the arbitrary and severe approximations, doubts re- 
main about their conclusiveness. 

Our works starts from recognizing that chemical inter- 



actions are ruled by various selection rules, and for the 
present situation spin-conservation M is expected to be 
relevant. In gas-phase chemistry it is well known that O2, 
when in its triplet ground state, is rather inert when the 
other reactant and the product are spin singlets. Inter- 
estingly, this role of the O2 spin has not attracted much 
attention in the 02/Al(lll) interaction, although it was 
e.g. studied for the adsorption of oxygen on Si(lOO) 
The appropriate theoretical modeling should then con- 
strain the spin to the O2 Hilbert subspace, preventing 
charge transfer, as well as spin quenching before the sys- 
tems interact. Such a spin-constrained DFT approach 
has neither been formulated nor evaluated for molecule- 
surface scattering so far. We will show that it not only 
gives a good description of the sticking coefficient (cf. 
Fig. ^ empty circles), but may also explain the enig- 
matic abstraction mechanism. 

Let us briefly describe the theoretical method enabling 
us to study the dynamics of an O2 molecule that re- 
mains in its spin-triplet configuration. Only very close 
to the surface transitions to other configurations of the 
02/Al(lll) system may set in. In order to calculate the 
spin-triplet PES we follow the work of Dederichs et al. 
jlSj l , for which one must first define the Hilbert subspace 
of the O2 molecule. As the DMoP code employs an atom- 
centered basis set, we use for this all basis functions that 
are also needed to provide a good description of the free 
O2 molecule. Then, for any position of the O2 molecule, 
we request that the total electron spin in this Hilbert sub- 
space is one. In practice this approach involves the self- 
consistent filling of the four partial densities of states of 
the spin-up and spin-down, O2 and Al(lll) sub-systems. 
This is formulated in terms of an auxiliary field in order 
to properly include the effect of the spin-constraint on 
the total energy 0. 

Before discussing the results obtained with this ap- 
proach, we remind of two general problems of present-day 
Kohn-Sham-DFT: First, even with gradient corrected xc 
functionals the description of the binding energy of the 
free O2 molecule is rather bad. Going from the O2 spin- 
triplet ground state to two free oxygen atoms, each of 
them also in the spin-triplet ground state, the errors of 
our calculated binding energies with respect to the ex- 
perimental value (5.1 eV 0) are: 2.3 eV (LDA), 1.0 
eV (PBE), 0.6 eV (BLYP), and 0.5 eV (RPBE). Fortu- 
nately, for the part of the PES, that is important for 
the sticking coefficient, we find that different function- 
als give results that differ by much less, indicating some 
favorable error cancellation. Below we will therefore re- 
strict our discussion to the PBE and the RPBE, since 
they represent the extreme cases for the gradient cor- 
rected functionals, yielding the strongest and smallest 
overbinding in the O2 molecule, respectively. A second 
noteworthy problem arises because the expectation value 
of is not defined in Kohn-Sham-DFT. For the present 
case this implies that the multiplet structure is not well 



3 



z(A) 
(a) 5.0 




z(A) 
(b)5.o 




z(A) 

(C)5.0 



4.0 



3.0 



2.0 



1.0 



Y lop view 2 side view 




i.oJ ^ ■ ^ 1.0^ 

1.0 1.4 1.8 2.2 r (A) 1.0 1.4 1.8 2.2 r (A) 1.0 1.4 1.8 2.2 r (A) 

z (A) z (A) z (A) 

(d)5.0| III h|' I'M' ; '. ' ' I (e)5.0| I II ' I 'I | [» ' ' ■ ' . ' ■ : ' . ■ (f)5.0T 

^ ^ ^^H,. ton view side view ^ ^ Ion view ^ side view ^ 





FIG. 2: Two-dimensional (elbow) cuts through the six-dimensional PESs calculated for three different situations, always using 
DFT-RPBE (see text): adiabatic (a,d), triplet (6, e) and singlet PES (c, /). The energies are shown as a function of the O2 
bond length r and of the distance Z of the O2 center of mass from the surface. The angles and lateral positions are indicated 
in the insets. The energy zero (green/yellow border) corresponds to a free triplet O2 molecule. Contour lines are drawn at 
0.2 eV intervals. Dissociation barriers (if present) are labeled (eV). 



described 0, 0. In free O2 the many-body ground 
state belongs to the triple degenerate state which 
is followed by two singlets, namely a doubly degenerate 
"'^Ag level (0.98 eV above the ground state), and a non 
degenerate ^Sg level (1.63 eV above the ground state). 
While the total energy of the spin-triplet ground state 
is described well, the ^Ag and ^Sg states are not de- 
scribed appropriately, since here DFT with jellium-based 
xc functionals describes a certain mixture of multiplets. 
A reasonable approximation to the true spin-singlet state 
is instead obtained by a spin-unpolarized calculation 0, 
which is for PBE 1.1 eV (for RPBE 1.2 eV) higher than 
the spin-triplet ground state. 

Figure [3 shows two cuts through the calculated six- 
dimensional PESs for three situations: the adiabatic ap- 
proximation (discussed in the introduction), the spin- 
triplet PES (using constrained DFT) and the spin- 
unpolarizcd calculation, which is the best we can do to 
describe the spin-singlet PES. Whereas the two elbow 
plots of the adiabatic PES (cf. Fig. [5] left panels) do 
not exhibit sizeable energy barriers toward dissociative 
adsorption, we find clear barriers on the triplet PES (cf. 
Fig. El middle panels). In fact, inspecting the whole six- 
dimensional triplet PES there is always an energy barrier 
(the lowest one is 0.05 eV). The right panels of Fig. [21 



show the corresponding cuts through the singlet PES, 
which never exhibits any energy barriers. Clearly, an O2 
molecule prepared in the singlet state would therefore 
react most efficiently with the Al(lll) surface. Since 
the spin forbidden transition to the triplet ground state 
can only proceed by scattering with another molecule, 
the long lifetime of a singlet O2 should render molecular 
beam experiments possible to verify this proposition. 

The sticking coefficient for these PESs is calculated as 
described above, i.e., using the "divide and conquer" ap- 
proach H, 13 ■ The results for the adiabatic and the 
triplet PESs, using the RPBE functional, are given in 
Fig. n Obviously, the spin-triplet PES gives a stick- 
ing curve in good agreement with the experimental re- 
sult. However, when the O2 and Al(lll) wave func- 
tions overlap at close distances, spin transfer will occur 
with a certain probability. Due to the uncertainty in 
the description of the singlet-PES, it is at present not 
very meaningful to perform a quantitative evaluation of 
these transition probabilities. A rough estimate of the 
importance of transitions bringing the system away from 
the triplet-PES is instead provided by the width of the 
2tt* Kohn-Sham resonance, which is the level that car- 
ries the spin. At large distance the width is zero, and 
it gradually increases upon approach to the surface. For 
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FIG. 3: Potential energy along the reaction path shown as 
dashed line in Figs. b, c (solid line = triplet PES, dotted 
line = singlet PES, dashed line = adiabatic PES). The red ar- 
row indicates the classical trajectory of a thermal O2 molecule 
constrained to the triplet PES, with CTP marking the classi- 
cal turning point. At this point the coupling, represented by 
the width of the O 2p Kohn-Sham level (dash-dotted line), is 
only just emerging. 



a one-dimensional cut through configuration space this 
is shown in Fig. 13 The peak width remains quite nar- 
row and even at the point where the triplet and singlet 
PESs cross it is only about 0.1 eV. In general, the life- 
time of the 27r* electrons should be compared to the time 
the molecule spends between the classical turning point 
(CTP) and ca. 5 A away from the surface. For ther- 
mal molecules (cf. the arrow and the CTP point in Fig. 
13) the comparison is: lifetime 3 ps vs. time of pres- 
ence « 1 ps. We therefore conclude that for thermal O2 
molecules (and even for all molecules with a kinetic en- 
ergy below 0.2 eV) transitions away from the triplet 
PES will not play a big role. Our results then suggest 
that particularly these lowest energy molecules should be 
repelled by the barriers on the triplet PES, well before 
there is significant hybridization of wave functions, i.e. 
before relaxation towards the adiabatic ground state oc- 
curs. Only for higher kinetic energies, transitions will 
gradually become important, leading to higher sticking 
coefficients than in the "theory-triplet" curve shown in 
Fig. n We also note that the PESs of the PBE and 
RPBE functionals are similar, but quantitatively differ- 
ences exist. These differences have noticeable influence 
on the calculated sticking curve only for kinetic transla- 
tional energies below 0.2 eV. As the RPBE gives a better 
description for free O2 we place a higher credibility on 
its PES. Details will be discussed elsewhere 0. 

Analyzing the approaching O2 molecule in greater 
detail reveals finally another interesting feature. For 
molecules that approach in an orientation perpendicu- 
lar to the surface (or close to this) the spin is shifted 
to the atom that is further away from the surface. We 
believe this to be the onset of adsorption by the abstrac- 



tion mechanism. In this way one O atom can adsorb in 
a singlet state, while the spin is efficiently carried away 
with the other O atom that is either repelled back into 
the vacuum or to a distant place at the surface. Calcu- 
lating the full dynamics of this process, i.e. going beyond 
the onset of dissociation important for the sticking coef- 
ficient, requires the explicit consideration of forces on the 
Al atoms, which we are implementing at present. 

In summary, we have shown that spin selection rules 
can play an important role for O2 scattering at metals. 
They imply that O2 molecules should travel in a spin- 
triplet configuration up to distances close to the surface 
where hybridization with metal-surface states becomes 
significant. This is particularly important for systems 
with a low DOS at the Fermi level; for transition metals 
we expect that the high density of c?-states at the Fermi 
level can more easily take up the spin. At Al(lll) spin 
selection leads to a very low sticking probability for ther- 
mal O2 molecules in the triplet ground state, while O2 
molecules prepared in the singlet configuration should 
adsorb with high probability. Similar effects as those 
discussed in this paper should just as well play a role for 
other substrates with a low jellium-like density of states 
at the Fermi level, and for other molecules with a high- 
spin ground state. 
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